Quỹ đạo là gì? Các nghiên cứu khoa học liên quan đến Quỹ đạo

Quỹ đạo là đường chuyển động của một vật thể quanh vật thể khác dưới tác dụng của lực hấp dẫn hoặc tương tác, có thể là tròn, elip hay hyperbol. Trong cơ học lượng tử, quỹ đạo mô tả vùng xác suất xuất hiện của hạt như electron quanh hạt nhân, gọi là orbital, mang tính phi định xứ và xác suất.

Khái niệm quỹ đạo

Quỹ đạo (orbit) là đường chuyển động mà một vật thể mô tả khi di chuyển dưới ảnh hưởng của một trường lực, phổ biến nhất là lực hấp dẫn. Trong cơ học cổ điển, quỹ đạo mô tả hành vi động học của các thiên thể như hành tinh, vệ tinh và sao chổi, thường có hình elip, tròn hoặc hyperbol tùy theo điều kiện động năng và lực tác động. Trong lĩnh vực lượng tử, khái niệm quỹ đạo cũng được sử dụng, nhưng mang nghĩa xác suất phân bố vị trí của các hạt vi mô như electron xung quanh hạt nhân nguyên tử.

Khái niệm quỹ đạo xuất hiện từ thời cổ đại trong thiên văn học cổ điển, nhưng chỉ được mô tả chính xác về mặt toán học nhờ công trình của Johannes Kepler và Isaac Newton vào thế kỷ 17. Quỹ đạo không chỉ mô tả chuyển động theo không gian ba chiều mà còn chịu ảnh hưởng bởi thời gian và năng lượng của hệ thống. Ở cấp độ vi mô, quỹ đạo không còn mang tính hình học chính xác mà trở thành khái niệm toán học trừu tượng, gắn liền với hàm sóng trong cơ học lượng tử.

Ngày nay, quỹ đạo là công cụ lý thuyết và thực nghiệm quan trọng trong các ngành như vật lý thiên văn, cơ học vũ trụ, hóa học lượng tử và công nghệ hàng không vũ trụ. Sự hiểu biết về quỹ đạo cho phép con người dự đoán chính xác chuyển động của vật thể trong vũ trụ, thiết kế lộ trình cho vệ tinh, tàu thăm dò và các thiết bị không gian khác.

Phân loại quỹ đạo theo vật lý học

Quỹ đạo có thể được phân loại theo hệ quy chiếu vật lý mà trong đó chúng được nghiên cứu hoặc ứng dụng. Việc phân loại này giúp xác định chính xác mô hình toán học và công cụ phân tích phù hợp cho từng lĩnh vực. Trong vật lý học, ba nhóm quỹ đạo chính được công nhận gồm quỹ đạo thiên văn, quỹ đạo nhân tạo và quỹ đạo nguyên tử.

  • Quỹ đạo thiên văn: Là quỹ đạo tự nhiên của các thiên thể quay quanh một vật thể khối lượng lớn hơn, như các hành tinh quay quanh Mặt Trời, mặt trăng quay quanh hành tinh chủ, hoặc sao chổi tiếp cận hệ Mặt Trời.
  • Quỹ đạo nhân tạo: Là các đường bay của vệ tinh nhân tạo hoặc tàu vũ trụ, được điều khiển bởi các cơ quan không gian. Quỹ đạo này có thể được lập trình và điều chỉnh theo thời gian.
  • Quỹ đạo nguyên tử (orbitals): Trong vật lý lượng tử, không phải là đường đi vật lý mà là không gian trong đó xác suất hiện diện của electron đạt cực đại. Các orbital có thể có hình cầu (s), hình số 8 (p), hoặc các dạng phức tạp hơn như d, f.

Phân loại quỹ đạo không chỉ giúp thuận tiện trong phân tích mà còn hỗ trợ thiết kế nhiệm vụ khoa học hoặc kỹ thuật phù hợp. Ví dụ, các vệ tinh viễn thông sử dụng quỹ đạo địa tĩnh để đảm bảo tín hiệu ổn định, trong khi vệ tinh quan sát Trái Đất cần quỹ đạo đồng bộ Mặt Trời để thu thập dữ liệu nhất quán theo thời gian.

Các định luật mô tả quỹ đạo

Quỹ đạo vật thể được mô tả chính xác thông qua các định luật cơ bản trong cơ học thiên thể. Hai hệ thống lý thuyết chủ yếu chi phối là định luật Kepler (mô tả hình học và thời gian của quỹ đạo) và định luật vạn vật hấp dẫn của Newton (mô tả lực chi phối quỹ đạo). Theo Kepler, mọi hành tinh chuyển động quanh Mặt Trời theo quỹ đạo elip với Mặt Trời nằm ở một tiêu điểm. Diện tích quét được bởi đường nối hành tinh với Mặt Trời là như nhau trong các khoảng thời gian bằng nhau.

Về lực hấp dẫn chi phối quỹ đạo, định luật Newton phát biểu rằng mọi cặp vật thể khối lượng m1 và m2 cách nhau một khoảng r sẽ hút nhau với lực: F=Gm1m2r2F = G \frac{m_1 m_2}{r^2} trong đó GG là hằng số hấp dẫn vũ trụ. Công thức này cho phép xác định chính xác gia tốc và biến thiên vận tốc của vật thể trong quỹ đạo.

Ở cấp độ vi mô, các quy tắc cổ điển không còn áp dụng được. Trong cơ học lượng tử, phương trình Schrödinger đóng vai trò xác định “quỹ đạo” lượng tử của các hạt như electron. Không tồn tại đường đi rõ ràng, thay vào đó là xác suất phân bố hàm sóng: H^Ψ=EΨ\hat{H} \Psi = E \Psi Trong đó H^\hat{H} là toán tử Hamilton mô tả năng lượng hệ thống, Ψ\Psi là hàm sóng, và EE là năng lượng riêng tương ứng.

Quỹ đạo tròn, elip và hyperbol

Hình dạng quỹ đạo của một vật thể phụ thuộc vào năng lượng toàn phần và mô men động lượng của nó. Dựa vào độ lệch tâm (ee) của quỹ đạo, người ta phân loại thành ba dạng chính: tròn, elip và hyperbol. Độ lệch tâm là đại lượng mô tả mức độ lệch khỏi hình tròn của một quỹ đạo conic.

  • Quỹ đạo tròn: Có độ lệch tâm e=0e = 0, là trường hợp lý tưởng hóa trong đó lực hướng tâm và tốc độ cân bằng tuyệt đối.
  • Quỹ đạo elip: 0<e<10 < e < 1, là dạng phổ biến nhất trong thiên văn học. Mọi hành tinh trong hệ Mặt Trời đều di chuyển theo quỹ đạo elip với mặt trời nằm tại một tiêu điểm.
  • Quỹ đạo hyperbol: e>1e > 1, xảy ra khi vật thể có động năng đủ lớn để thoát khỏi ảnh hưởng hấp dẫn của thiên thể trung tâm. Thường áp dụng cho sao chổi hoặc tàu vũ trụ bay qua hành tinh mà không quay lại.

Bảng sau tóm tắt các dạng quỹ đạo theo độ lệch tâm và ứng dụng thực tiễn:

Loại quỹ đạo Độ lệch tâm (e) Đặc điểm Ứng dụng/Hiện tượng
Tròn 0 Bán kính không đổi, lực hấp dẫn cân bằng tốc độ Vệ tinh viễn thông
Elip 0 < e < 1 Hai tiêu điểm, tốc độ biến thiên theo vị trí Hành tinh quay quanh sao
Hyperbol e > 1 Chuyển động mở, vật thể rời khỏi hệ Sao chổi, tàu thăm dò liên sao

Việc tính toán hình dạng quỹ đạo là yếu tố sống còn trong thiết kế lộ trình bay cho tàu vũ trụ. Các mô hình toán học và phần mềm như GMAT hoặc STK được sử dụng để mô phỏng các tham số quỹ đạo và tối ưu hóa tiêu hao năng lượng.

Quỹ đạo ổn định và bất ổn định

Quỹ đạo ổn định là dạng quỹ đạo mà vật thể quay quanh một trung tâm hấp dẫn trong thời gian dài mà không bị thay đổi đáng kể về hình dạng hoặc vị trí. Ổn định này thường xảy ra trong các hệ hai vật có tỷ lệ khối lượng chênh lệch lớn như Trái Đất và Mặt Trăng, hoặc Mặt Trời và các hành tinh. Trong các hệ như vậy, lực hấp dẫn trung tâm gần như chi phối hoàn toàn quỹ đạo, và không có yếu tố gây nhiễu đáng kể.

Tuy nhiên, không phải tất cả quỹ đạo đều ổn định. Một quỹ đạo có thể trở nên bất ổn định khi bị ảnh hưởng bởi các yếu tố ngoài như nhiễu loạn hấp dẫn từ vật thể thứ ba, lực kéo khí quyển, áp suất bức xạ mặt trời hoặc các nhiễu loạn do sai số trong điều khiển. Trong các mô hình nhiều vật thể (three-body or n-body problem), việc xác định sự ổn định của quỹ đạo trở nên cực kỳ phức tạp và thường phải dựa vào mô phỏng số hoặc lý thuyết hỗn loạn (chaotic dynamics).

Bảng sau phân biệt một số đặc điểm giữa quỹ đạo ổn định và bất ổn định:

Tiêu chí Quỹ đạo ổn định Quỹ đạo bất ổn định
Thời gian tồn tại Rất dài (hàng triệu năm) Ngắn hoặc trung bình
Ảnh hưởng từ môi trường Không đáng kể Cao
Khả năng dự đoán Cao Khó dự đoán (chaotic)
Ứng dụng Quỹ đạo hành tinh, vệ tinh địa tĩnh Quỹ đạo bay qua, chuyển tiếp Hohmann

Ứng dụng của quỹ đạo trong kỹ thuật và khoa học

Hiểu rõ về quỹ đạo là điều kiện tiên quyết trong các lĩnh vực như hàng không vũ trụ, viễn thông, định vị toàn cầu và khoa học vũ trụ. Các vệ tinh viễn thông được đặt trên quỹ đạo địa tĩnh (GEO) để giữ vị trí cố định so với một điểm trên Trái Đất, cho phép truyền tín hiệu ổn định và liên tục. Trong khi đó, hệ thống định vị toàn cầu (GPS) sử dụng nhiều vệ tinh trên quỹ đạo trung bình (MEO) đồng bộ để cung cấp dữ liệu thời gian và vị trí cho người dùng toàn cầu.

Quỹ đạo cũng được áp dụng trong các sứ mệnh liên hành tinh như tàu vũ trụ sử dụng quỹ đạo chuyển tiếp Hohmann để tiết kiệm nhiên liệu khi di chuyển giữa hai quỹ đạo khác nhau. Trong thiết kế tên lửa và tàu vũ trụ, việc tính toán đúng quỹ đạo không chỉ tiết kiệm năng lượng mà còn đảm bảo an toàn cho thiết bị và con người.

  • Vệ tinh viễn thông: Đặt trên quỹ đạo GEO với chu kỳ 24 giờ, độ cao ~35.786 km.
  • GPS: Sử dụng 24–32 vệ tinh MEO ở độ cao khoảng 20.200 km.
  • Sứ mệnh không gian: Thiết kế theo các cửa sổ phóng tối ưu hóa năng lượng, áp dụng quỹ đạo elip hoặc hyperbol.

Thông tin chi tiết về mạng lưới vệ tinh TDRS của NASA có thể xem tại NASA TDRS.

Quỹ đạo trong vật lý lượng tử

Trong cơ học lượng tử, khái niệm “quỹ đạo” không còn là một đường đi vật lý xác định. Electron trong nguyên tử không di chuyển theo đường tròn hay elip như hành tinh, mà được mô tả bởi “orbital” – vùng không gian nơi có xác suất tìm thấy electron là cao nhất. Orbital không chỉ có kích thước mà còn có hình dạng đặc trưng, phụ thuộc vào các số lượng tử n (chính), l (xung lượng), m (từ) và s (spin).

Hàm sóng Ψ\Psi của electron được tính từ phương trình Schrödinger, và xác suất hiện diện là: P(x)=Ψ(x)2P(x) = |\Psi(x)|^2 Các orbital như s (hình cầu), p (hình số 8), d và f có cấu trúc phức tạp, cho phép giải thích các hiện tượng hóa học như liên kết, lai hóa và phổ hấp thụ. Do ảnh hưởng của nguyên lý bất định Heisenberg, ta không thể đồng thời xác định chính xác vị trí và động lượng của electron.

So với quỹ đạo cổ điển, quỹ đạo lượng tử mang tính xác suất và phi định xứ, là nền tảng cho toàn bộ hóa học hiện đại, đặc biệt trong tính toán cấu trúc phân tử và thiết kế vật liệu điện tử.

Quỹ đạo trong lý thuyết tương đối

Thuyết tương đối rộng của Einstein đã thay đổi hoàn toàn cách nhìn về quỹ đạo trong môi trường hấp dẫn mạnh. Theo thuyết này, vật thể không di chuyển theo lực hấp dẫn như trong Newton, mà theo các đường trắc địa (geodesic) trong không–thời gian cong. Lực hấp dẫn không còn là lực thực sự mà là biểu hiện của độ cong không–thời gian.

Quỹ đạo trong không–thời gian cong có thể giải thích được hiện tượng mà cơ học Newton không lý giải được, ví dụ như dịch chuyển cận điểm quỹ đạo của Sao Thủy, hiệu ứng thấu kính hấp dẫn, hoặc sự uốn cong ánh sáng gần hố đen. Một ví dụ quan trọng là việc đo được thời gian chậm đi khi ở gần vật thể có khối lượng lớn – một yếu tố cần hiệu chỉnh trong hệ thống GPS hiện đại.

Tham khảo bài viết gốc của Viện Vật lý Hoa Kỳ tại AIP - Origins of General Relativity.

Phân tích quỹ đạo bằng công cụ mô phỏng

Với sự phát triển của công nghệ tính toán, các công cụ mô phỏng quỹ đạo đã trở thành công cụ thiết yếu trong thiết kế, điều khiển và dự đoán chuyển động của các thiết bị không gian. Các phần mềm như GMAT (General Mission Analysis Tool) của NASA, STK (Systems Tool Kit) của AGI, hoặc MATLAB Aerospace Toolbox cho phép người dùng mô phỏng lực, gia tốc, biến đổi quỹ đạo, và tối ưu hóa năng lượng cho các sứ mệnh không gian phức tạp.

Các yếu tố như nhiễu loạn hấp dẫn, lực cản khí quyển, độ lệch khối lượng hành tinh và áp suất bức xạ được đưa vào mô hình. Điều này giúp tính toán chính xác quỹ đạo thực tế, xác định thời gian cửa sổ phóng và điểm rơi của vệ tinh hoặc tên lửa.

  • GMAT: Mô phỏng quỹ đạo 3D, tối ưu hóa delta-v, phân tích hiệu suất tên lửa.
  • STK: Dùng trong mô hình hóa hệ thống vệ tinh, kết nối dữ liệu, và giám sát sứ mệnh.
  • SPICE Toolkit: Cung cấp tập dữ liệu về hành tinh và quỹ đạo được NASA sử dụng trong các sứ mệnh không gian.

Tài liệu tham khảo

  1. Roy, A. E. (2005). Orbital Motion (4th ed.). CRC Press.
  2. Vallado, D. A. (2013). Fundamentals of Astrodynamics and Applications. Microcosm Press.
  3. Murray, C. D., & Dermott, S. F. (1999). Solar System Dynamics. Cambridge University Press.
  4. Griffiths, D. J. (2018). Introduction to Quantum Mechanics (3rd ed.). Cambridge University Press.
  5. Misner, C. W., Thorne, K. S., & Wheeler, J. A. (1973). Gravitation. W.H. Freeman.
  6. NASA. (2024). TDRS: Tracking and Data Relay Satellite System. Link
  7. AIP. (2023). Origins of General Relativity. Link

Các bài báo, nghiên cứu, công bố khoa học về chủ đề quỹ đạo:

Các phương pháp quỹ đạo phân tử tự nhất quán. XX. Một tập hợp cơ sở cho hàm sóng tương quan Dịch bởi AI
Journal of Chemical Physics - Tập 72 Số 1 - Trang 650-654 - 1980
Một tập hợp cơ sở Gaussian loại thu gọn (6-311G**) đã được phát triển bằng cách tối ưu hóa các số mũ và hệ số ở cấp độ bậc hai của lý thuyết Mo/ller–Plesset (MP) cho trạng thái cơ bản của các nguyên tố hàng đầu tiên. Tập hợp này có sự tách ba trong các vỏ valence s và p cùng với một bộ các hàm phân cực chưa thu gọn đơn lẻ trên mỗi nguyên tố. Tập cơ sở được kiểm tra bằng cách tính toán cấu ...... hiện toàn bộ
#cơ sở Gaussian thu gọn #tối ưu hóa số mũ #hệ số #phương pháp Mo/ller–Plesset #trạng thái cơ bản #nguyên tố hàng đầu tiên #hàm phân cực #lý thuyết MP #cấu trúc #năng lượng #phân tử đơn giản #thực nghiệm
Phương pháp quỹ đạo phân tử tự nhất quán. XII. Phát triển bổ sung bộ cơ sở dạng Gaussian cho nghiên cứu quỹ đạo phân tử của các hợp chất hữu cơ Dịch bởi AI
Journal of Chemical Physics - Tập 56 Số 5 - Trang 2257-2261 - 1972
Hai bộ cơ sở mở rộng (được gọi là 5–31G và 6–31G) bao gồm các hàm sóng nguyên tử được biểu diễn dưới dạng kết hợp tuyến tính cố định của các hàm Gaussian được trình bày cho các nguyên tố hàng đầu từ cacbon đến flo. Những hàm cơ sở này tương tự như bộ 4–31G [J. Chem. Phys. 54, 724 (1971)] ở chỗ mỗi lớp vỏ hóa trị được chia thành các phần bên trong và ngoài được mô tả tương ứng bằng ba và mộ...... hiện toàn bộ
#quỹ đạo phân tử #hàm cơ sở Gaussian #cacbon #flo #năng lượng tổng #cân bằng hình học #phân tử đa nguyên tử
Phương pháp quỹ đạo phân tử tự trùng khớp: Mở rộng cơ sở kiểu Gaussian cho nghiên cứu quỹ đạo phân tử của các phân tử hữu cơ Dịch bởi AI
Journal of Chemical Physics - Tập 54 Số 2 - Trang 724-728 - 1971
Một tập hợp cơ sở mở rộng của các hàm số nguyên tử được biểu diễn dưới dạng các tổ hợp tuyến tính cố định của các hàm Gaussian được trình bày cho hydro và các nguyên tố hàng đầu tiên từ cacbon đến flo. Trong tập này, được mô tả là 4–31 G, mỗi lớp vỏ bên trong được đại diện bởi một hàm cơ sở duy nhất được lấy từ tổng của bốn hàm Gaussian và mỗi quỹ đạo hoá trị được tách thành các phần bên t...... hiện toàn bộ
#Hàm Gaussian #cơ sở dữ liệu phân tử #ổn định cấu trúc #tối ưu hóa năng lượng #quỹ đạo phân tử
Quản trị hợp tác trong lý thuyết và thực tiễn Dịch bởi AI
Journal of Public Administration Research and Theory - Tập 18 Số 4 - Trang 543-571 - 2008
Tóm tắt Trong vài thập kỷ qua, một hình thức quản trị mới đã xuất hiện để thay thế các cách làm chính sách và thực thi theo kiểu đối kháng và quản lý. Quản trị hợp tác, như đã được biết đến, mang lại sự tương tác giữa các bên công và tư qua các diễn đàn tập thể cùng với các cơ quan công để tham gia vào việc ra quyết định định hướng theo đồng thuận. Trong bài báo nà...... hiện toàn bộ
#quản trị hợp tác #mô hình tùy biến #ra quyết định đồng thuận #xung đột #hợp tác #niềm tin #lãnh đạo #thiết kế thể chế
Phân tích quỹ đạo liên kết tự nhiên của dimers nước gần-Hartree-Fock Dịch bởi AI
Journal of Chemical Physics - Tập 78 Số 6 - Trang 4066-4073 - 1983
Chúng tôi đã tiến hành phân tích quỹ đạo liên kết tự nhiên của liên kết hydro trong dimers nước đối với hàm sóng gần-Hartree–Fock của Popkie, Kistenmacher, và Clementi, mở rộng các nghiên cứu trước đây dựa trên các tập hợp cơ sở nhỏ hơn và hình học ít thực tế hơn. Chúng tôi phát hiện rằng các tương tác có thể được mô tả một cách chính xác là ''chuyển giao điện tích'' (đặc biệt là tương tác...... hiện toàn bộ
#liên kết hydro #quỹ đạo liên kết tự nhiên #dimers nước #tương tác điện tích #van der Waals
Liệu pháp quay điều biến thể tích: IMRT trong một quỹ đạo gantry đơn Dịch bởi AI
Medical Physics - Tập 35 Số 1 - Trang 310-317 - 2008
Trong nghiên cứu này, một nền tảng tối ưu hóa kế hoạch điều trị mới được trình bày, nơi liệu pháp được thực hiện một cách hiệu quả và chính xác trong một quỹ đạo quay điều chế động duy nhất. Những cải tiến trong chăm sóc bệnh nhân đạt được thông qua việc định vị hình ảnh và điều chỉnh kế hoạch đã dẫn đến sự gia tăng thời gian điều trị tổng thể. Liệu pháp bức xạ điều biến cường độ (IMRT) cũ...... hiện toàn bộ
Các con đường phát thải toàn cầu dưới các kịch bản kinh tế xã hội khác nhau để sử dụng trong CMIP6: một tập dữ liệu về các quỹ đạo phát thải đã được hài hòa đến cuối thế kỷ Dịch bởi AI
Geoscientific Model Development - Tập 12 Số 4 - Trang 1443-1475
Tóm tắt. Chúng tôi trình bày một bộ chín kịch bản về các quỹ đạo phát thải trong tương lai từ các nguồn nhân tạo, một sản phẩm chính của thí nghiệm ScenarioMIP trong CMIP6. Kết quả của các mô hình đánh giá tích hợp cho 14 loài phát thải khác nhau và 13 lĩnh vực phát thải được cung cấp cho mỗi kịch bản với những chuyển tiếp nhất quán từ dữ liệu lịch sử được sử dụng trong CMIP6 sang các quỹ ...... hiện toàn bộ
Nhận thức Đạo đức trong Các Tổ Chức Kinh Doanh: Ảnh Hưởng của Các Yếu Tố Liên Quan đến Vấn Đề và Bối Cảnh Xã Hội Dịch bởi AI
SAGE Publications - Tập 53 Số 7 - Trang 981-1018 - 2000
Sự nhận thức của cá nhân về các vấn đề đạo đức là bước đầu tiên quan trọng trong quá trình ra quyết định đạo đức. Dựa vào nghiên cứu trong lĩnh vực nhận thức xã hội và đạo đức kinh doanh, chúng tôi đã giả thuyết rằng sự nhận thức đạo đức bị ảnh hưởng bởi các yếu tố liên quan đến vấn đề (mức độ hậu quả của vấn đề đạo đức và cách trình bày vấn đề theo cách đạo đức) và các yếu tố liên quan đ...... hiện toàn bộ
#Nhận thức đạo đức #yếu tố liên quan đến vấn đề #bối cảnh xã hội #ra quyết định đạo đức #nghiên cứu thực địa
Sự sửa đổi quỹ đạo của chuyển động chỉ trỏ để đáp ứng sự thay đổi vị trí mục tiêu Dịch bởi AI
Journal of Neurophysiology - Tập 49 Số 2 - Trang 548-564 - 1983
1. Các chuyển động tay hướng về một mục tiêu liên quan đến chuyển động tại khớp vai và khớp khuỷu tay và bị hạn chế trong mặt phẳng sagittal đã được nghiên cứu. Trong một số chuyển động, vị trí mục tiêu thay đổi đột ngột, do đó yêu cầu một sự điều chỉnh có chủ ý của quỹ đạo bởi đối tượng con người. 2. Thời gian phản ứng để điều chỉnh quỹ đạo tương đương với thời gian phản ứng để khởi động ...... hiện toàn bộ
Hiệu chuẩn hình học trước phóng và trên quỹ đạo của VIIRS Suomi NPP Dịch bởi AI
Journal of Geophysical Research D: Atmospheres - Tập 118 Số 20 - 2013
Tóm tắtCảm biến VIIRS (Visible Infrared Imaging Radiometer Suite) được phóng vào ngày 28 tháng 10 năm 2011 trên vệ tinh Suomi National Polar‐orbiting Partnership (SNPP). VIIRS có 22 dải quang phổ covering from 0.412 µm đến 12.01 µm, bao gồm 16 dải độ phân giải trung bình (dải M) với độ phân giải không gian danh nghĩa 750 m tại điểm thẳng đứng, năm dải độ phân giải ...... hiện toàn bộ
Tổng số: 463   
  • 1
  • 2
  • 3
  • 4
  • 5
  • 6
  • 10